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Abstract 

Functionalized polcyclic aromatic hydrocarbons (PAHs) are an interesting class of molecules in which the 
electronic state of the graphene-like hydrocarbon part is tuned by the functional group. Searching for new 
types of donor and acceptor molecules, a set of new PAHs has recently been investigated experimentally 
using ultraviolet photoelectron spectroscopy (UPS). In this work, the electronic structure of the PAHs is 
studied numerically with the help of B3LYP hybrid density functional. Using the ASCF method, electron 
binding energies have been determined which affirm, specify and complement the UPS data. Symmetry 
properties of molecular- orbitals are analyzed for a categorization and an estimate of the related signal 
strength. While cr-like orbitals are difficult to detect in UPS spectra of condensed film, calculation provides 
a detailed insight into the hidden parts of the electronic structure of donor and acceptor molecules. In 
addition, a diffuse basis set (6-31 1++G**) was used to calculate electron affinity and LUMO eigenvalues. 
The calculated electron affinity (EA) provides a classification of the donor/acceptor properties of the studied 
molecules. Coronene-hexaone shows a high EA, comparable to TCNQ, which is a well-known classical 
acceptor. Calculated HOMO-LUMO gaps using the related eigenvalues have a good agreement with the 
experimental lowest excitation energies. TD-DFT also accurately predicts the measured optical gap. 

Keywords: PAHs, charge transfer complex, electron binding energy, vertical ionization, Koopmans' theorem, UPS, 
electron affinity 



* |naghavi@uni-mainz.de 



2 



I. INTRODUCTION 



Charge transfer salts have become of particular interest during the past two decades, in which 
both the interest in electronics on the molecular scale as well as the methods of chemical synthesis 
has strongly increased. Typically, charge transfer (CT) complexes, consisting of yr-electron donors 
(D) and acceptors (A), show one or several properties like metallic conductivity, superconductiv- 
ity, and magnetism that are relevant for nanoelectronic applications yj. 



conductors and superconductors are 
TTF) donor and various acceptors 12, 



20. Many of the organic 
)ased on the bis(ethylenedithio)tetrathiafulvalene (BEDT- 
,|3|]. Furthermore, poly cyano derivatives have mainly been 
used as acceptors like tetracyanoquinodimethane (TCNQ) [141]. Recently, a novel type of donor and 
acceptor molecules based on poly cyclic aromatic hydrocarbons (PAHs) has been investigated ||5|]. 
In this article, a detailed numerical study of the electronic properties of new PAHs is presented, 
which allows to identify the experimental results and complements the picture of experimentally 
determined orbital states. 

Flat aromatic molecules as polycyclic aromatic hydrocarbons (PAHs), have attracted attention 
because of their extended 7r-electron systems |^ 7]. Garito and Heeger concluded that the exten- 
sion of the TT system would lead to a lowering of intramolecular Coulomb repulsion in the anions of 
the acceptors, resulting in more stable radical anions and in highly conducting CT complexes [SI]. 
Larger PAH molecules called nanographene have been synthesized in the bottom-up approach 



with desired size, structure, and symmetry The electronic structure of PAHs can be tuned via 
the type and number of substitutions in periphery. This way the functionalized PAHs may serve 
as a prototype for the creation of a new class of molecules with tailored chemical and electronic 
properties. However, although the electronic structure of different donor and acceptor molecules 
have been widely studied, little attention has been paid by theory and experiment to illustrate the 
electronic structure and the effect of functional groups on donors and acceptors based on PAHs. 

In this study we calculate the electronic structure of new acceptors and donors based on 
coronene and pyrene (see Fig. [T]). A hybrid density functional (B3LYP) method is used to cal- 
culate the ionization potentials (IPs). The results are compared with ultraviolet photoelectron 
spectroscopy (UPS) measurements [|5D. The comparison allows to assign the measured signals to 
specific molecular orbitals, providing a self consistent validation of the experimental results. Fur- 
thermore, the calculations provide the IPs of orbitals that could not be measured, experimentally. 

Unfortunately, gas phase measurements were not possible in the experiment because of the 



small amount of material available (a few milligrams). Instead, UPS measurements had been 
performed on PAH multilayers on a Au surface. As the molecules are preferentially planar to the 
surface while measurements are made in the perpendicular direction a-type orbitals could not be 
measured in practice. In the numerical study all orbitals can be investigated so that the respective 
IPs become accessible. In general, the IPs of PAHs adsorbed on a metal surface are shifted by 
a constant energy value with respect to the gas phase results. The discrepancy stems from the 
photo-hole screening induced by the image charge in the metal Isj]. A similar effect occurs in 
multilayer films in which the screening is induced by the neighboring molecules. The energy shift 
between the theoretical results that corresponds to isolated molecules or the dilute gas phase and 
the measurements on the condensed multilayers can be taken into account by adjusting the HOMO 
levels of both systems. Applying this energy shift to all IPs must provide a good correspondence 
between the experimentally measurable orbitals and the theoretical results, which serves as a use- 
ful consistency check of both measurement and calculation. For the given set of molecules the 
energy shift is very small so that the IP results for isolated molecules and those in the multilayer 
can be compared directly. 

Coronene has a high symmetry (Dg/i), a small band gap and the right size for processing tech- 
niques. Therefore, coronene and its derivatives are promising for applications like molecular 
electronics optoelectronics or sensors [Il0l - ll3n . Kato et al. [114,] also predicted superconductiv- 
ity of coronene in a charge transfer configuration. In this article, properties of hexamethoxy- 
coronene (HMC) and coronene-hexaone (CHO) are studied. The functional groups are tailored 
such that HMC is a donor while CHO is an acceptor molecule. Furthermore, calculations of 
tetramethoxypyrene (TMP) and pyrene-tetraone (PTO) have been performed and compared with 
experimental results. These poly aromatic molecules have a different symmetry and structure in 
comparison with coronene based molecules, which makes a comparison of these two families 
instructive. Calculation of tetracyanoquinodimethane (TCNQ), a well-known classical acceptor, 
is used as a benchmark for the calculations of the acceptor molecules. All investigated molecules 
are shown in Figure [IJ Electronic structures of the donor and acceptor molecules are theoretically 
studied by B3LYP/6-31G(d) level of theory and compared with the ultraviolet photoemission 
spectroscopy (UPS) measurements that were carried out in parallel [|5D. Electron binding energies 
or ionization potentials of molecules are calculated with the delta self consistent field (ASCF) 
method using the B3LYP hybrid functional. 
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FIG. 1 . Structure of the donor and acceptor molecules studied in this work. 



n. COMPUTATIONAL DETAILS 



The calculations were performed using the GaussianOS code 11 1511 a nd visualized by VMD lll6n 



Calculations were carried out with the B3LYP hybrid functional 
tionals combine the Becke's three parameter exchange potentia' 



tionals of Lee- Yang-Parr 



17h. This combination of func- 



19|1 with the correlation func- 



21D and Vosko-Wilk-Nusair [|22D. B3LYP is one of the often em- 



ployed hybrid functional [|230 used in theoretical studies of molecules. 

The standard 6-31G(d) basis set ll24ll is employed for geometry optimizations and energy cal- 
culations. This level of theory ensures a reliable estimate of structural, energetic and spectro- 



scopic parameters for PAHs, which requires an extremely high computational effort for the large 
molecules. 

Electron binding energies or ionization potentials of molecules are calculated by the delta self- 
consistent field (ASCF) method applied together with the B3LYP hybrid functional. Vertical ion- 
ization potentials are calculated using ASCF method in order to compare with UPS data in which 
vertical ionization potentials (IP^) are measured. In IP^, final and ground state have the same ge- 
ometry. First vertical ionization potentials (IP^,i) are calculated as the difference in total energy 
between a molecule and its radical cation at the former geometry. For the higher ionization ener- 
gies (IPt,,i+„), the energy difference (A E) between IP„ i and the HOMO energy is added to the 



energy of next higher orbitals Shomo-u as a uniform shift [|25l - l27ll 



IPv,i - \eHOMo\ = A^ (1) 

l^HOMO-nl + = IPy,i+n (2) 

According to Eq. [2] the eigenvalue of the HOMO-n orbital should be used in order to calculate 
the electron binding energy (BE) of HOMO-n. A standard view expressed in the literature for 
Kohn-Sham (KS) orbital energies of the occupied molecular orbitals lower than HOMO, is that 
they are merely auxiliary quantities. However, several authors have pointed to the interpretative 



power of KS orbitals in traditional qualitative molecular orbital schemes [|27 



also show a uniform shift with respect to ionization energies [|27l 



M, and KS orbitals 



311]. In this sense, B3LYP with 



its exact Hartree-Fock exchange contribution leads to better description of ionization energies by 



reducing the related shift [SOO . 

The B3LYP/6-31 l-i-i-G(d,p) level of theory was implemented on both the anion and the neutral 
species for the electron affinity calculations. In case of an anionic final state, the usage of dif- 
fuse function is crucial to ensure the accuracy of the calculation. The spin unrestricted B3LYP 
(UB3LYP) is employed in which species are ionized with an open-shell doublet (one unpaired 
electron) electron configuration. Spin contamination is found in accepted limits for the open shell 
systems, with <S^> values of about 0.75-0.78 in all cases. At the stationary points a vibrational 
frequency calculation is performed for each numerically optimized structure to confirm that it is 
located at the minimum of the potential energy surface. 

The results obtained with the ASCF method are compared with Hartree-Fock calculations 



based on the traditional Koopmans' theorem approach [|32h . According to Koopmans' theorem, 



the ionization potential (IP) is the sign-reversed orbital energy of Hartree-Fock eigenvalues: IP 



-Ehomo- Although KT is based on some simplification, it shows a reasonable agreement with ex- 
perimental data in several cases 113 311 . Koopmans' theorem neglects the relaxation effect by using 
the frozen-orbital approximation. However, this error is frequently compensated by the oppositely 
directed error due to the electron correlation effect, neglected in the Hartree-Fock (HF) method. 



Therefore, the Koopmans' theorem is a crude but useful and fast approach [13311 . 

Moreover, eigenvalues based on B3LYP/6-311-i-+G(d,p) are used to estimate HOMO-LUMO 
gaps and LUMO energies. HOMO-LUMO gap based on the energy difference between the HOMO 
and LUMO eigenvalues are compared with the measured optical gaps. TD-B3LYP/6-31G(d) is 
also used to estimate the optical band gap (referring to the singlet excited state). As shown by 
Zhang et al the first excitation energy from Tp-B3LYP calculations leads to a good prediction of 
the optical gap for most of the DFT methods [34]. Moreover, HOMO-LUMO gaps obtained with 
TD-DFT are not very sensitive to the basis set and it is shown that 6-31G(d) can show reasonable 

n 

agreement for optical properties with the experimental results for the hydrocarbon molecules B4I - 



3611 . The energy difference between the HOMO and LUMO eigenvalues and TD-DFT show the 



similar results and both are in an excellent agreement with the experimental results. 



III. RESULTS AND DISCUSSION 



The result and discussion section is divided into six parts. The first part is about the experi- 
mental setup in UPS. In the second part, the symmetry of molecules and consequences of orbital 
symmetries on the UPS experiment are described. The effect of functional groups on coronene 
and pyrene derivatives is described in the third part. In the following two parts, experimental and 
theoretical results of the binding energies for coronene and pyrene derivatives are presented. In 
the sixth part, the calculated LUMO and electron affinity (EA) values for the studied compounds 
are discussed. 



A. UPS measurements 

Since only small amounts of the new molecules were avaible, UPS measurements have not been 
performed on molecules in the gas phase but on UHV-deposited molecules forming some multi- 
layers on a Au surface. In the measurement, the work function of a clean Au surface increases 
significantly if the acceptor TCNQ or coronene-hexaone is adsorbed (isjl. It can be understood in 
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terms of a partial charge transfer from the metal surface to the first molecular layer thus forming 
the interface dipole with its negative side pointing toward the organic film. Opposite behavior is 
observed for the donors HMC, TMP and coronene. The decrease of the work function of about 
0.5 eV indicates the Pauli push-back effect leading to a push-back of the split-out electrons and 
thus to a reduction of the surface dipole of the metal. In practice, due to charge transfer between 
molecules and surface, the polar metal-organic interfaces change the work function and thus the 
binding energies with respect to vacuum level. 

The binding energies in condensed films of the studied acceptors and donors (Fig.[T]) are mea- 
sured by UPS [|5D. Li the UPS spectra the signals derived from the occupied frontier orbitals of 
the organic overlayer are referenced to the Fermi energy of the Au substrate. By adding the value 
for the work function of the molecular film, the binding energies are referenced to the vacuum 
level. An acceptor molecule in direct contact with a metal surface may form a negatively charged 
species as was found for TCNQ on Ni(l 1 1) ll37n . TCNQ on Au, on the other hand, was found to be 
almost neutral jsS] because its acceptor character does not suffice to form a strong charge-transfer 
complex with Au, which is the metal with the highest Pauling electronegativity. For coronene on 
a Au surface an ionization potential of IP=7.1 eV is measured (see Table I]) which is very close to 
the reported value of 7.29 eV for coronene in the gas phase 113911 . As a result, one has that the Au 
surface shows neutrality with respect to both donor (coronene) and acceptor (TCNQ) molecules 
and experimental and theoretical IP values can be compard without considering a large energy 
shift. 

The following work function values have been measured for films with a thickness of several 
layers deposited on Au: coronene-hexaone (5.6 eV), hexamethoxycoronene (4.8 eV), coronene 
(4.8 eV), tetramethoxypyrene (4.9 eV) and TCNQ (5.9 eV). 

Electron binding energies for the molecular orbitals (MOs) can be extracted from the UPS spec- 
tra by assigning peak positions of the spectra, but one by one assignment for all UPS signals is not 
possible since in some cases they result from superposition of several occupied molecular orbitals 
with closely spaced binding energies. However, in most cases the HOMO signal is well separated 
from the deeper-lying levels and can be evaluated with high precision. 
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B. Symmetry analysis 



For the analysis and comparison of UPS data with the calculated ionization potentials, the sym- 
metry of the molecular orbitals must be extracted, a or tt symmetry is determined, followed by 
related degeneracy (E, A, ...).This information is shown in Figure[2]for the frontier orbitals and in 
Table IJ and In] for the lower lying orbitals of studied molecules. 

Having n symmetry is essential in order to observe the orbital in UPS. The oxygen-derived 
a states appear only weakly in photoelectron spectroscopy because of the lower cross section of 



y, 



0-2p in comparison with carbon tt at the photon energy (21.2 eV) used in the experiment 
Another reason for the weak a-derived signals in the experiment might be that the molecules are 
predominantly oriented parallel to the surface and the UPS spectrometer probes the normal emis- 
sion direction. Electronic structure calculations allow to determine the electron binding energies 
of the a-type orbitals, which are not accessible in the experiment. This way a detailed picture of 
the electronic structures of donor and acceptor molecules is available. Orbital symmetries with a 
and TT character of MOs are listed in Tables Uln] and Figure |2l 

Coronene and its derivatives have a degeneracy of two for the HOMO and all orbitals with E 
symmetry because of more than a 2-fold axis (coronene has 6-fold and HMC and CHO have 3-fold 
axis). CHO and HMC have a doubly degenerated HOMO with a tt characteristic as depicted in 
Figure[2l For HMC and TMP there is no symmetry plane, so a and it orbitals are not clearly distin- 
guishable for the lower lying orbitals. However, they are reported according to major contribution 
of cr- or TT-type symmetry in the orbital's notations in Tables H] and HIl 

Pyrene derivatives and TCNQ have a non-degenerate orbital because of the 2-fold axis. Pyrene- 
tetraone (PTO) has an exceptional situation because HOMO and HOMO- 1 have cr-type symmetry 
and are in-plane as shown in Figure [B Therefore, they are not traceable in UPS because of the low 
cross section of 0-2p since these MOs are located mainly on the oxygen atoms (see Figure[2l). For 
PTO, reliable UPS binding energies could not be determined due to charging. 

To sum up, a survey of Tables |I] and |II] reveals that the presence of O atoms and cyanid groups 
increases the a-type part of the occupied orbitals of acceptor molecules. Thus, acceptors have 
more cr-symmetry orbitals near the HOMO. The experimentally inaccessible a-type orbitals are 
provided by the electron structure calculations as shown in Tables HI HIl 



Coronene (Dsh) 




HOMO (E2u) LUMO (Eig) 

Hexamethoxycoronene (HMC) (D3) 




Coronene-hexaone (CHO) (Dsh) Tetracyanoquinodimethane (TCNQ) (D2h) 




HOMO(E") LUMO(Ag) H0M0(B3„) LUMO (Bag) 



Tetramethoxypyrene (TMP) (D2) Pyrene-tetraone (PTO) (Dah) 




HOMO(Bi) LUMO (A) HOMO (Bzu) LUMO (Bg^) 



FIG. 2. HOMO and LUMO orbitals computed for the studied molecules. The symmetry of each molecule 
and the symmetry of the related molecular orbitals are shown in parentheses. Among all PAHs donor and 
acceptor molecules just PTO has a cr-type HOMO orbital. All the coronene derived molecules show a 
doubly degenerate HOMO. 
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C. Effect of functional groups 



The electron-withdrawing carbonyl groups of CHO and PTO reduce the charge density in 
the TT system, while in HMC and TMP, electrons are donated to the tt system by the methoxy 
(OMe) groups. Therefore, the aromatic rings in CHO and PTO are positively charged, while the 
methoxy substituted molecules are negatively charged. Thus, CHO is an acceptor in comparison 
with coronene and HMC, and PTO is an acceptor compared to TMP. 

The bond lengths of CHO along the symmetry plane are distinctly larger than those of coronene. 
The first C-C bond length at the edge of the double bonded oxygen atoms increases from 1 .37 A 
(coronene) to 1.54 A for CHO. The other bond lengths increase by about 0.02 A with respect to 
coronene. For CHO, all the bond lengths inside the aromatic ring are larger than those in HMC and 
coronene. In HMC, the C-C bond lengths are of the same size as those of coronene, but MuUiken 
population shows that the C atoms of the aromatic rings are slightly more charged compared to 
coronene. In principle, the similar behavior is found for TMP and PTO. 

Donor and acceptor molecules are best investigated by analyzing the molecular orbitals. Any 
change in the functional groups of the molecules will be reflected in the binding energies and the 
molecular orbital energies. Therefore, evaluation of molecular orbitals shows how the functional 
groups tune the electronic structure. Also, with the help of ab-initio calculations one can analyze 
details of the electronic structure and in particular a orbitals which are not accessible by UPS 
measurement. Also underlying orbital structure of overlapping UPS signals can be analyzed by 
the calculation. 

D. Electron binding energies of coronene derivatives 

Theoretical and experimental binding energies of coronene and its derivatives coronene- 
hexaone (CHO) and hexamethoxycoronene (HMC) are listed in Table HI A survey of Table |I] 
reveals that vertical electron binding energies (IPs) calculated in ASCF-B3U^P approximation 
are, on the whole, in good agreement with the corresponding experimental IPs. Calculations 
based on the ASCF method correspond to the vertical IP of the molecules in the gas phase. In 
general, they may differ from measurements of condensed layers on a metal surface by a constant 
energy value. As shown in Table H] the energy shift between the multilayer measurements and the 
calculations corresponding to the gas phase are very small, which corresponds to the fact that the 
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IP=7.1 eV of coronene on Au is close to the reported gas phase value IP=7.29 eV [!39i]. 

For coronene and the derivatives the Koopmans' energies for the HOMO differ by less than 
0.5 eV from the ASCF binding energies. For the deeper-lying orbitals, the Koopmans' energies 
are too large in comparison with ASCF and experimental results. This deviation is caused by 
correlation and relaxation effects which are neglected in Koopmans' theorem. For comparison, 
results from Koopmans' theorem are presented for HMC in Table HI For coronene and CHO devi- 
ations of the Koopmans' energies are of the same order of magnitude. Note that the cr-like orbitals 
of CHO appear only very weakly in UPS spectra. ASCF-B3LYP results can cover a wide range 
of MOs and show good agreement with the UPS data. 

The first ionization potential (binding energy of HOMO) of HMC is shifted up by 0.5 eV in 

TABLE I. Electron binding energies (eV) for coronene and its derivatives with the symmetry labels cal- 
culated by B3LYP/6-31G* and compared to vertical excitation energies from ultra violet photoelectron 
spectroscopy (UPS) and Koopmans' theorem (KT). 



CH0(D3/,) HMC(D3) Coronene(D6/») 



MOs 


Symm. 


ASCF 


UPS 


Symm. 


KT 


ASCF 


UPS 


Symm. 


ASCF 


UPS 


HOMO 


E"(vr) 


8.78 


8.8 


E(^) 


6.86 


6.41 


6.5 


E2M(vr) 


6.84 


7.1 


HOMO-1 


E'(a) 


8.80 




E(7r) 


8.74 


7.66 




Eig(7r) 


8.11 




HOMO-2 


A[{a) 


8.82 




A2(7r) 


9.25 


7.80 


>7.9 


B25(vr) 


8.56 


>%.% 


HOMO-3 


A'/(vr) 


9.23 




Ai(7r) 


9.44 


8.15 




Bi3(7r) 


8.62 




HOMO-4 


E'(CT) 


10.05 




E((T) 


11.03 


8.38 




A2„(vr) 


9.84 




HOMO-5 


A'2(CT) 


10.11 




Ai(cT) 


11.59 


8.51 






9.94 




HOMO-6 


E"(7r) 


10.33 




E(^) 


11.65 


9.11 




E2«(vr) 


10.25 








>10.3 












HOMO-7 


A^'(^) 


10.70 




A2(Tr) 


11.65 


9.17 


>9.3 


B2«(ct) 


11.09 




HOMO-8 


E"(7r) 


11.80 




A2(7r) 


12.47 


9.52 




El^^(c^) 


11.13 




HOMO-9 


E'(a) 


11.84 




E(7r) 


12.64 


9.71 




A2g(cT) 


11.54 




HOMO- 10 


A'/(vr) 


12.37 




E(7r) 


12.94 


9.95 




B1^.(CT) 


11.55 





comparison to coronene. As displayed in Table H hexamethoxycoronene has three main signals in 
the experimental spectrum below 10 eV [sl. The first of them at 6.5 eV is the HOMO with vr-like 
symmetry with the double degenerate representation E (Figure[2l). The calculation yields a binding 
energy (A-SCF) of 6.41 eV in excellent agreement with the photoelectron spectra. The calculated 
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binding energy of electrons in the HOMO of HMC is about 0.4 eV lower than that of coronene 
due to the six methoxy groups on the periphery which increase the charge density in the it system. 
The effect of the methoxy groups is weaker than that of electron-withdrawing keto groups of 
CHO since the methoxy groups contain the electron-donor methyl and the acceptor oxygen which 
partially compensate each other. There are three pairs of methoxy groups in the ring system and 
each pair is in trans configuration (see Figure [H and [2l). The methoxy groups can easily rotate in 
a solvent but this rotation will be blocked on the surface. The orientation of the methoxy group 
configurations has a strong influence also on the interaction between metal surface and molecule. 
The trans configuration of methoxy groups may reduce the contact of the aromatic ring system 
with the surface so that it can effectively suppress charge transfer between molecule and metal 
surface. The second signal in the UPS spectra appears at 7.9 eV and is assigned to the group 
HOMO-1, HOMO-2 and HOMO-3 all having vr-like symmetry. Note that there is no strict vr- and 
(T-type symmetry in HMC, but we can identify V-like' and 'a-like' symmetry by comparison with 
the parent molecule coronene. HOMO-4 and HOMO-5 are very weak in the experimental spectra 
because they have ci-like symmetry, and contain a high contribution from oxygen 2p which has a 



41|1 . Another reason 



low cross section at the photon energy (21.2 eV) used in the experiment [ 4o[ 
for the weak a-derived signals in the experiment might be that the molecules are predominantly 
oriented parallel to the surface and the UPS spectrometer probes the normal emission direction. 
For perfectly parallel orientation, a emission is forbidden by symmetry selection rules. The UPS 
signal at 9.3 eV corresponds to the group of the HOMO-6 to HOMO- 8 states with vr-like symme- 
try. 

Coronene-hexaone (CHO) constitutes a more complicated adsorption behavior than HMC. 
CHO is a planar molecule with high electron affinity. The first layer being in direct contact with 
surface experiences a substantial charge transfer from Au showing up in terms of a new interface 
state i^. For comparison with theory we refer the UPS data for higher CHO coverage, where the 
interface state has disappeared in the spectra. In the multilayer coverage two signals rise in the 
region below 1 1 eV. As shown in Table HI the first signal corresponds to the HOMO and appears 
at 8.8 eV in very good agreement with the ASCF calculation (8.78 eV). This value is similar to 
TCNQ which is a strong classical acceptor(Figure The binding energy of electrons in the 
HOMO of CHO has increased by more than 2 eV with respect to coronene because of the elec- 
tronegative keto groups. HOMO-1 and HOMO-2 as well as HOMO-4 and HOMO-5 are weak 
signals in the experimental spectra since they are derived from the oxygen 2px,y orbitals that are 
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aligned in the molecular plane. HOMO-6 and HOMO-7 (tt symmetry) constitute the second UPS 
signal at 10.3 eV in the measured spectra. The ASCF calculation again shows good agreement. 
Obviously the six carbonyl groups strongly alter the electronic structure of CHO compared to 
unfunctionalized coronene, and make CHO a strong acceptor. 

E. Binding energies of pyrene derivatives 

For tetramethoxypyrene (TMP) and tetracyanoquinodimethane (TCNQ), the theoretical 
binding energies are smaller than the experimental values. For the HOMO of TMP the difference 
is 0.42 eV, similar as for the overlapping HOMO-1 and HOMO-2 signals. HOMO-3 to HOMO-4 
and HOMO-9 to HOMO- 10 have a-like symmetry and appears weakly in the experimental spectra 
due to a large oxygen 2p content. The 7r-like group HOMO-6 to HOMO-8 can not be resolved 
experimentally and shows up as one signal centered at 9.2 eV. 

TABLE II. Same as table U but for pyrene derivatives and TCNQ. 



PT0(D2J TMP(D2) TCNQ(D2/J 



MOs 


Symm. 


ASCF UPS 


Symm. 


ASCF 


UPS 


Symm. 


ASCF 


UPS 


HOMO 


B2«(fT) 


8.70 


Bi(^) 


6.18 


6.7 


B3«(^) 


8.87 


9.3 


HOMO-1 




8.90 


BsCvr) 


7.09 




Big(7r) 


10.49 












>1.1 




aO.8 


HOMO-2 




8.92 


B2(^) 


7.30 




B29(vr) 


10.91 




HOMO-3 


A«(7r) 


9.50 


A(^) 


8.03 




B39(fT) 


11.56 




HOMO-4 


Big(^) 


9.67 


B3(CT) 


8.06 






11.58 




HOMO-5 




9.96 


A(7r) 


8.28 




A„(cr) 


11.87 




HOMO-6 


Bl„((j) 


10.26 


Bi(7r) 


8.66 




Bi9(7r) 


11.88 




HOMO-7 


B3^.(vr) 


10.60 


B2(vr) 


8.72 


>9.3 


Bi„(7r) 


11.89 




HOMO-8 


B25(vr) 


11.66 


B2(^) 


8.96 






12.01 




HOMO-9 


BagCcT) 


12.08 


Bi(^) 


9.42 




B39(fT) 


12.34 




HOMO- 10 


Ag(a) 


12.10 


B2(CT) 


9.71 




B2«(<t) 


12.38 





For pyrene-tetraone (PTO) reliable UPS binding energies could not be determined due to 
charging. Moreover, the HOMO has a symmetry and is thus not observable in the UPS experiment. 
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Charging was also observed for thick CHO films (above 10 nm thickness). For TCNQ, HOMO- 
1 and HOMO-2 agree fairly well with the experimental signal centered at 10.6 eV. HOMO-3 to 
HOMO-5 have a symmetry and appear only weakly in UPS spectra. 



F. LUMO and electron affinity 

The strength of an acceptor molecule is measured by its electron affinity (EA) which is the 
energy released when adding one electron to the lowest unoccupied molecular orbital. An accep- 
tor must have a high electron affinity, while low electron affinity is required for effective donor 
molecules. Since EA depends on the lowest unoccupied molecular orbital (LUMO), it can be cal- 
culated more accurately by considering the final state (anion). In order to describe the anion with 
an extra electron sufficiently, a diffuse basis set is required. We use the B3LYP/6-31 1++G(d,p) 
basis to study donor and acceptor molecules. 

The calculated EA of coronene is 0.48 eV which is in good agreement with the reported ex- 



perimental value of 0.470±0.009eV ll42ll . A high electron affinity is the major characteristic of 



strong acceptors like TCNQ. Reported values for the EA of TCNQ vary with the applied mea- 



surement technique [|43l-l46l1 . Calculations values of EA are in the range of 3.4±4eV [47]. Using 
B3LYP/6-31 1-i-i-G**, we obtained an EA of 3.62 eV, while B. Milian et al. report a value of 
3.42 eV calculated with B3LYP/cc-pVDZ and 3.28 eV calculated with CCSD/aug-cc-pVDZ. All 



iQH. 



these values are close to the reported experimental value of EA=3.3±3 11441 

The electron affinity of CHO is 3.5 eV (see first column of Table HIH) which is close to the 
electron affinity of TCNQ (3.62 eV). Therefore CHO is expected to be a strong electron acceptor. 
CHO has a large electron affinity and an extended tt system which makes it an interesting alterna- 
tive to common acceptors. 

In Table Hill the experimental gaps are compared with the theoretical gaps. To estimate the 
optical gap, both TD-DFT and the difference of B3LYP HOMO-LUMO eigenvalues are used and 
the results are listed in the second and third column of Table Hill respectively. In the case of large 
molecules there is a good agreement between the eigenvalue differences and the measured optical 
gap. In the earlier study, it is shown that there is a large deviation between the calculated gap 



using B3LYP and experimental results [13311 . In the same paper, discrepancy depends to the size of 
molecule and it decreases by increasing the size of hydrocarbons. Therefore, one can concludes 
the good agreement in our calculations should be attributed to the large size of the hydrocarbons. 
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TABLE III. ASCF electron affinity (EA) using 6-311++G** basis set. Moreover the HOMO-LUMO 
gaps are calculated using TD-DFT/6-31G* (singlet excitation), second column and B3LYP/6-31 1++G** 
(LUMO-HOMO eigenvalue differences) third column. All values are given in eV. 



Molecule 


EA 


EA 


Gap 


Gap 


Gap 




(B3LYP) 


(Exp) 


(TD-DFT) 


(£lumo — £homo) 


(Exp) 


Coronene 


0.48 


0.47^ 


3.99 


3.23 


3.54^ 


Coronene-Hexaone 


3.50 




2.98 


2.57 


3.2(5^ 


Hexamethoxycoronene 


0.51 




3.81 


3.09 


3.30^ 


Pyrene-Tetraone 


2.48 




3.48 


2.22 




Tetramethoxypyrene 


0.41 




3.67 


3.58 




TCNQ 


3.62 


3.3^ 


3.04 


2.54 


3.10^ 



''Ref 1 


42I1 




''Ref 1 




46] 


'^Ref 1 






•iRef 


m 





For all compounds, the calculated HOMO-LUMO eigenvalue differences are smaller than the mea- 
sured gaps. The gaps obtained from TD-DFT are spread symmetrically around the experimental 
values. In general, the calculated gaps from TD-DFT and the HOMO-LUMO eigenvalues deviate 
less than 20% from the experimental gaps, which is a reasonable result for this quantity. The fair 
agreement between the experimental values and the results calculated with B3LYP/6-31G* indi- 
cates that this calculation method may be used also for the larger molecules of this family called 
nanographene for which more correlated methods like post Hartree-Fock or even DFT with the 
diffuse basis sets are difficult or not applicable. 



IV. SUMMARY 

In this paper, we have reported theoretical and experimental electron binding energies of dif- 
ferent MOs for coronene, pyrene and their methoxy and keton-derivatives as polycyclic aromatic 
hydrocarbons and compared them with TCNQ as a well-known classical acceptor. The studied 
compounds demonstrate how the functional groups can tune the molecules' electronic state, pro- 
viding moderate donor or strong acceptor properties. 

Comparison of electron structure calculations and photoelectron spectroscopy data provides 
new insights into the electronic structure of the new compounds. The calculated molecular or- 
bitals have been used to identify and complement the experimental results. UPS measurements 
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of the compounds on a Au surface provide IPs in an agreement with experimental and calculated 
values for the gas phase. The absence of an energy shift, induced by the metal substrate or the 
neighboring molecules, is so far not understood, but the good agreement between the various UPS 
signals and the calculated vr-type orbitals demonstrates the direct comparability of the molecule's 
IPs in the condensed and the gas state. In fact, ASCF-B3LYP reproduces the whole range of IPs 
in the measured UPS spectra, while Koopmans' theorem shows a good agreement just for the first 
IP. The calculated electron affinities allow to rate the donor or acceptor character of the molecules. 
They confirm that CHO is a strong acceptor. Its high electron affinity (3.5 eV), which characterizes 
a good acceptor, is almost as large as the affinity of TCNQ (3.62 eV). HMC, coronene and TMP 
with their low ionization potentials and rich aromatic systems are suitable candidates for donors. 
HMC is an appropriate choice for surface applications due to its trans isomerism. Especially inter- 
esting for the formation of charge transfer complexes are the molecules CHO and HMC. Our re- 
sults show that they are good acceptor and donor molecules, respectively, and with their coronene 
centers they have a very similar molecular structure. This is analogous to the TTF-TCNQ pair 
with the benzene core which is studied since a long time. The proposed coronene-based donor- 
acceptor molecules could be the starting point of a new generation of charge transfer salts based 
on molecular nanographenes. 
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